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ABSTRACT 
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pLn ' We present high resolution 890/im images of two 20 mJy submilhmeter galax- 

ies, SMMJ123711+622212 and MIPS J142824.0+352619, obtained using the Sub- 
millimeter Array (SMA). Using submillimeter interferometric observations with 
an angular resolution of 2'.'5, the coordinates of these high redshift sources are de- 
^ ■ termined with an accuracy of 0'.'2. The new SMA data on SMMJ123711+622212 

^ ■ reveal an unresolved submm source offset to the east by O'/S from an optical galaxy 

^ ■ found in deep HST images, suggesting either a large galaxy with a dusty cen- 

^ ' tral region, or an interacting galaxy system. The SMA image of hyper-luminous 

O ■ (Lfir = 3.2 X 10^3 L0) source MIPS J142824.0+352619 provides a firm upper 

• limit to the source size of < 1'.'2. This constraint provides evidence that the 

foreground lens is only weakly affecting the observed high FIR luminosity. 

O 

^ ■ Subject headings: galaxies: formation, galaxies: starburst, cosmology: observa- 

tions, galaxies: high redshift, submillimeter 
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1. Introduction 

The discovery of high redshift submilhmeter sources has significantly improved our un- 
derstanding of the star formation history in the early universe. Negative k-correction allows 
observation of the thermal dust emission at 850/im almost independent of the redshift up 
to 2; ~ 10 (Blain et al. 2002). Deep observations using the SCUBA bolometer on the James 
Clerk Maxwell Telescope (JCMT) have unveiled the presence of distant submm sources 
(Small, Ivison & Blain 1997; Barger et al. 1998; Hughes et al. 1998; Eales et al. 2000; Cowie, 
Barger & Kneib 2002; Scott et al. 2002; Borys et al. 2003; Webb et al. 2003; Wang, Cowie & 
Barger 2004). The primary origin of the submm emission is believed to be the reprocessed 
dust emission from newborn stars in young galaxies. While these discoveries are attended by 
a great number of high resolution follow-up optical/NIR imaging studies, the 14" resolution 
of the JCMT at 850/xm yields a large error circle which is too coarse for a precise deter- 
mination of the optical/NIR counterparts to these sources. Deep optical imaging typically 
shows several optical/NIR sources within the SCUBA beam. To date, the most successful 
ways to obtain precise astrometry on the target are to obtain high resolution, deep 1.4 GHz 
radio images (Ivison et al. 1998; Barger, Cowie Sz Richards 2000; Chapman et al. 2001; 
Ivison et al. 2002; Dunlop et al. 2004), or to obtain interferometric 1.3 mm continuum im- 
ages (e.g. Downes et al. 1999). The former, however, does not identify robust optical/NIR 
counterparts for all of the sources, revealing counterparts for ~ 75% of the S'gso/im > 5 mJy 
sources with 5'i.4ghz > 30 yuJy (e.g. Ivison et al. 2002; Chapman et al. 2003; Borys et al. 
2004; Greve et al. 2004; Wang, Cowie & Barger 2004). Precise astrometry obtained using 
mm or submm interferometers can unambiguously identify the correct counterpart for the 
remaining radio-faint sources, and sources with multiple radio counterparts. High angular 
resolution submm observations also allow us to understand the true nature of the submm 
sources with established optical counterparts in which gravitational lensing is a possibility. 

We present recent Submilhmeter Array (SMA; Ho, Moran & Lo 2004) detections of 
two S85oiim ~ 20 mJy sources, SMMJ12371 1+622212 (hereafter GN 20; Pope et al. 2005) 
and MIPS J142824.0+352619 (hereafter MIPS-J1428; Borys et al. 2005). The 20.3 mJy 
source GN 20 was discovered in the recent SCUBA observations of the GOODS North Field 
(Giavalisco et al. 2004). GN 20, with its lOo" 850;um detection (Pope et al. 2005) and 5(t 
1.3 mm detection at IRAM PdBI (Pope et al. in preparation) make this one of the strongest 
submm sources discovered to date. It has very weak radio and undetectable 450;um emission, 
suggesting it lies at high redshift. MIPS- J 1428 was discovered in the Spitzer MIPS images 
of the NDWFS Bootes field (Jannuzi & Dey 1999; Soifer et al. 2004). It was detected 
(Borys et al. 2005) at 350//m using SHARC-II (Dowell et al. 2003) on the CSO, 1.4 GHz 
radio continuum at the VLA, and subsequent foUowup observations (Borys et al. 2005) 
revealed that it is an extremely luminous (3.2 ± 0.7) x 10^^ Lq starburst dominated galaxy 
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ai z — 1.325. Recent Keck-DEIMOS spectroscopy has revealed a, z — 1.034 galaxy directly 
aligned with MIPS-J1428 (Borys et al. 2005), possibly lensing it. High angular resolution 
submm observations allow us to determine whether or not the submm emission is coincident 
with the optical emission, and thereby search for any evidence for amplification of the FIR 
luminosity by gravitational lensing. We adopt Hq = 70 km s~^ Mpc"^, ilm = 0.3, fl^ = 0.7. 

2. Observation and Data Reduction 

GN 20 was observed on February 20 (track 1) and March 5, 2005 (track 2), and MIPS- 
J1428 was observed on March 8 (track 3) and April 4 (track 4), 2005 using 5-7 antennas in 
the compact configuration of the SMA. The SIS receivers were tuned to a center frequency 
of 342.883 GHz in the upper sideband (USB), yielding 332.639 GHz in the lower sideband 
(LSB). This tuning frequency was chosen to facilitate the receiver tuning and to achieve 
the optimal receiver performance. The target coordinates were obtained from the 1.3 mm 
IRAM PdBI detection of GN 20 {a (J2000) = 12'^37'"11.88^ S (J2000) = 62°22'12'.'00; 
Pope et al. in preparation), and the position of the proposed MIR galaxy of MIPS-J1428 
{a (J2000) = 14'^28"^24.10^ 5 (J2000) = 35°26'19'/00; Borys et al. 2005). All tracks were 
taken under good atmospheric opacity (i.e. T225 = 0.04 - 0.08). 

The SMA data were calibrated using the Caltech software package MIR, modified for 
the SMA. Antenna based passband cahbration was done using all of the planets and bright 
QSOs observed in a given track. For GN 20, antenna based time-dependent phase calibration 
was done using 1153+495, a 0.7 Jy QSO 14° away from the target. In addition, a 0.6 Jy 
QSO near GN 20, 1048+717 (14° from GN 20, 23° from 1153+495) was observed for a total 
of 12 minutes during each track. The detection of 1048+717 at the phase center empirically 
verifies and constrains the accuracy of the phase calibration referenced to 1153+495 (see §3). 
Similarly for the MIPS-J1428 tracks, two QSOs 1310+323 (0.6 Jy; 16° away from MIPS- 
J1428) and 1635+381 (1.0 Jy; 25° away from MIPS-J1428) were used together to calibrate 
the time dependent phase, and 1419+543 (0.4 Jy; 3° away from MIPS-J1428) was used to 
check the astrometry. Finally, absolute flux calibration was performed using Callisto and 
Mars. Imaging was carried out in MIRIAD (Sault, Teuben & Wright 1995). Maximum 
sensitivity was achieved by adopting natural weighting, which gave a synthesized beam size 
of 2'.'9 X 2'.'2 (P.A.= 12.2°) for GN 20 and 2'.'6 x 2'.'4 (P.A. = -47.2°) for MIPS-J1428. The 
rms noise after combining the two sidebands in two tracks was 2.1 mJy (GN 20) and 2.2 mJy 
(MIPS-J1428). 
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3. Results and Discussion 

3.1. SMMJ123711+622212 (GN 20) 

The SMA image of GN 20 is shown in Figure 1 (a). The real part of the visibility 
amplitudes do not decline as a function of projected baseline length, indicating that the 
emission is not spatially resolved (Figure 2). The upper limit on the source size is l'/2. The 
derived total flux from a point source model is 22.9 ±2.8 mJy, consistent with the 890;um flux 
of 18.1 mJy extrapolated from the SCUBA 850/im flux of 20.3 ±2.1 mJy (Pope et al. 2005). 
The derived coordinates for GN 20 are a (J2000) = 12'*37'"11.92", 5 (J2000) = 62°22'12'.'10, 
with statistical uncertainties in the flt of O'.'l for both a and 5. These arc consistent with 
estimated errors of ~ O'.'l 1 in a and ~ O'.'IS in 5 from a lOo" detection and a beam of 2'.'9 x 2'.'2. 

In order to check the robustness of our phase calibration and to estimate the systematic 
uncertainties in the SMA astromctry, we have imaged the test QSO (1048±717) using the 
same phase calibration we used to map GN 20. The resultant QSO map after adding the 
two sidebands and two tracks is shown in Figure 1 (a) (inset). A point source fit to the 
visibilities of 1048+717 gave a positional offset from the phase center of Aa = O'.'Ol ± 0'/02 
and = 0'.'06±0'.'02. The coordinates of this and all of the QSOs used in these observations 
were adopted from the Radio Reference Frame (Johnston et al. 1995), which are accurate 
to better than 3 mas. Hence the precise detection of 1048+717 at the phase center ensures 
that our phase calibration referenced to 1154+379 is robust. 

As an additional check, we fit a point source model to the visibihties of 1048+717 in 
each sideband of each track separately. The results show that the offsets from the phase 
center are consistent among the different sidebands and tracks in R.A. (Aa — O'.'Ol - 0'.'04), 
and slightly larger in Dec. for track 2 (A5 = 0'.'02 - 0'.'05), and a factor 10 larger in Dec. for 
track 1 (A5 ~ 0'.'30). The overall average offset from the phase center is Act = 0'.'02 ± 0.02 
and A5 = 0'.'16±0.14. The uncertainties in Dec. in track 2 are larger than the uncertainties 
in R.A. due to beam elongation in the north-south direction, while the source of the factor 
10 larger error in track 1 is not obvious from the data. These errors are consistent with 
the uncertainties of O'.'l - 0'.'15 expected from a maximum baseline error of 0.1 A at 230 GHz 
using a phase cahbrator that is 15° away from the target. These systematic tests prove the 
robustness of our phase calibration in each track and lend high confidence to the resulting 
positional accuracy of < O'.'l in R.A. and O'.'l - 0'.'2 in Dec. for GN 20. 

The new astrometric coordinates allow us to compare the submm source with high res- 
olution images in the publicly available deep Spitzer IRAC (Fazio et al. 2004) (Figure 1 (b)) 
and HST ACS (Ford et al. 1998) (Figure 1 (c)) images of the GOODS North Field. The 
absolute astrometric accuracy of both of these images is ~ O'.'l, and they are both tied to the 
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coordinate frame defined by the VLA positions in Richards (2000)^. The IRAC 3.6//m image 
reveals a source centered < 0'.'5 west of the submm coordinates, while the higher resolution 
HST V-haiid image reveals a faint optical source C'S to the west. From the analysis of avail- 
able ACS images, it is found that this optical source is a B-dropout galaxy {B = 27.2 ± 0.4, 
V = 25.2 ±0.1, i = 24A ± 0.1) which gives constraints on the probable redshift to be 
z ~ 3-4. We believe that the O'.'S offset between the SMA position and the HST position is 
significant, and there are several possible astrophysical explanations for this difference. The 
submm emission may arise from part of a large galaxy where y-band emission is completely 
obscured. Alternatively, GN 20 might be an interacting system where the optical galaxy is 
a companion to the dusty, more actively star forming galaxy. Although the observed FIR 
luminosities are significantly different, the apparent separation of these galaxies (~ 6 kpc 
at z ~ 3) suggest a close similarity to the Antennae system (NGC 4038/39), where the 
optical galaxies are separated by 7.5 kpc and most of the starburst activity is occurring in 
the medium between the two galaxies (e.g. Wang et al. 2004). 

GN 20 was suggested to be a two component source, GN 20.1 (20.3 mJy) and GN 20.2 
(11.7 mJy), separated by 18" (140 kpc) in the low resolution SCUBA image (Pope et al. 2005). 
The 890/im flux agreement between our SMA observation and the SCUBA measurement of 
GN 20.1 implies that GN 20 may be a two component source. However, the suggested 
position of GN 20.2 is beyond the half power point of the SMA primary beam, where the 
sensitivity is reduced by more than a factor of two. Detailed discussion in the context of 
multi-wavelength observations of this source will be provided in a forthcoming paper (Pope 
et al. in preparation). 



3.2. MIPS J142824.0-I-352619 (MIPS-J1428) 

The SMA map of MIPS-J1428 is shown in Figure 1 (d). The derived total fiux is 
18.4 ±2.5 mJy and shows excellent agreement with the 890/xm fiux of 19.5 mJy extrapolated 
fi-om the SCUBA 850/xm fiux of 21.9 ±1.3 mJy. MIPS-J1428 is not spatially resolved 
with the 2'.'5 beam (see Figure 2), and the derived coordinates from a point source fit are 
a (2000) = 14'*28'"24.06^ S (2000) = 35°26'19'.79, with uncertainties in the fit of O'.'l for 
both a and 6. The strong 890/xm detection allows us to make a higher resolution image 
using the visibilities of the longest baselines. The resulting 3a unresolved image provides a 
firm upper limit to the source size of < 1'.'2 (10 kpc aX z = 1.325), which is smaller than the 



^After applying the known positional offset in Dec. of C'SS between the GOODS-N images and the VLA 
catalog positions. See http://data.spitzer.caltech.edu/popular/goods/Documents/ for details. 
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size constraint given by the ~ 1'.'5 VLA 1.4 GHz resolution. 

The map of the test QSO 1419+543 after adding the two sidebands and two tracks is 
shown in Figure 1 (d) (inset). A point source fit to the visibihties gave a positional offset from 
the phase center of Aa = (y.'ll±0.05 and AS = (y.'23±0.05 for 1419+543. As with 1048+717 
for GN 20, we fit a point source model to the visibilities of 1419+543 in each sideband of 
each track separately. The results show that the offsets from the phase center had a wide 
range of values (Aa = 0'.'04 — 0'.'28 and A5 = C.'IS — 0'.'25) with the overall average offset 
from the phase center of Aa = 0'.'14 ±0.11 and A6 = 0'.'21 ± 0.05. These uncertainties are 
larger than those found in 1048+717 for GN 20, and arc slightly larger than the uncertainties 
expected from baseline errors. It is possible that other factors such as image smearing due 
to large phase noise may have introduced a small error in the source positions. Thus, we 
assess an uncertainty of ~ 0^'15 in R.A. and ~ 0'.'2 in Dec. for MIPS-J1428. 

Figure 1 (e) and (f) show the SMA 890/xm contours of MIPS-J1428 overlaid on the 
NOAO Deep Wide Field Survey K-band and I-band images. The astrometry of both of 
these images is tied to the reference frame defined by the USNO-A2.0 catalog, and the typ- 
ical rms of the residuals is 0'.'35 (Jannuzi et al. 2005). The optical/NIR galaxy seen here is 
unambiguously the galaxy aligned with the strong submm emission. The accurate astrome- 
try of the submm emission provided by the SMA allows tight constraints on the separation 
between the bright optical/NIR position, and rules out fainter IRAC sources detected nearby 
(Borys et al. 2005). Hence MIPS-J1428 could lie directly behind the foreground z = 1.034 
galaxy, which would potentially result in a large amplification of the submm source. Borys 
et al. (2005) use size/luminosity relationships to argue that despite the alignment, the am- 
plification is likely modest since the Einstein ring is of comparable size to the known physical 
scales of local ULIRGs. Assuming that the lensing is modest, the size is comparable to that 
of other high-redshift submm sources (Chapman et al. 2004). Using the star- formation rate 
estimated in Borys et al. (2005), our limit on the angular size, 9, of the object, we derive 
a lower limit on the SFR density of > 180(^/1. 2)|Li^^M0yr~^kpc~^ where n is the lensing 
amplification. This is larger than local ULIRGs (Mcurcr et al. 1997) and comparable to 
other high-redshift submm galaxies (Ghapman et al. 2004). 



4. Summciry 

We present SMA observations of two 20 mJy submm sources, GN 20 and MIPS-J1428. 
The positions of the submm sources are determined with O'.'l — 0'.'2 accuracy with these data, 
allowing precise identification of the correct optical galaxy counterpart to the bright submm 
emission in GN 20, and providing evidence that the foreground lens is only weakly affecting 
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the observed high FIR luminosity in MIPS-J1428. If many of the bright submm galaxies are 
slightly lensed objects similar to MIPS-J1428, then the implied source counts of inherently 
bright submm galaxies are over-predicted. Detailed studies, however, exist for only a few 
sources, and future surveys such as SHADES (Mortier et al. 2005) will provide important 
information about the star formation properties at high-redshifts. 

The Submillimctcr Array is a joint project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy and Astrophysics, and is funded 
by the Smithsonian Institution and the Academia Sinica. This work made use of observa- 
tions made with the Spitzcr Space Telescope, which is operated by JPL, California Institute 
of Technology under NASA contract 1407. This work made use of images provided by the 
NOAO Deep Wide-Field Survey (Jannuzi & Dey 1999) which is supported by the National 
Optical Astronomy Observatory (NOAO). NOAO is operated by AURA, Inc., under a co- 
operative agreement with the National Science Foundation. 
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Fig. 1.— (a) The SMA 890/im map of GN 20 and (mset) 1048+717. The lowest positive 
contours represent 2a, and the contours increase by la for GN 20, and by 2a for 1048+717. 
The lowest negative contour is 2a and increases by la. [b) The SMA 890/im contours 
overlaid on the Spitzer IRAC 3.6/im image and over the (c) HST ACS V^-band image, both 
obtained from the GOODS archive. The IRAC and HST images are corrected for the known 
0'.'38 offset in declination. The 4, 6, and 8 a contours from (a) are shown, and errorbars 
near the center show the astrometric accuracy of the SMA image, (d) The SMA map of 
MIPS-J1428 and (inset) 1419+543. The contours are the same as in GN 20. The smallness 
of the astrometric errors (O'.'l - 0'.'2) from the phase center in the 1048+717 and 1419+543 
maps prove the robustness of the astrometry of GN 20 and MIPS-J1428. The SMA 890/im 
map of MIPS-J1428 is shown overlaid on the NDWFS (e) K-band image and the (/) I-band 
image. The 3, 5, and 7 a contours from (d) are shown, and errorbars near the center show 
the astrometric accuracy of the SMA image. 
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Fig. 2. — The real part of the visibihty amphtudes vs. projected basehne length for (top) 
GN 20 and (bottom) MIPS-J1428. The solid horizontal lines represent a point source model 
with continuum fluxes of 22.9 mJy (GN 20) and 18.4 mJy (M1PS-J1428). The upper limit 
to the source sizes of both of these sources are comparable to the scale constrained by the 
longest baseline length which is ~ 1'.'2. 



